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arrays of coupled bistable elements. We present and compare experimental 
and numerical results for kink propagation under the influence of local and 



' global fluctuations. As demonstrated previously for local noise, an optimum 

' range of global noise power exists for which the medium acts as a reliable 

transmission "channel". We discuss implications for propagation failure in a 
model of cardiac tissue and present a general theoretical framework based on 
' discrete kink statistics. Valid for generic bistable chains, the theory captures 

o 

' the essential features observed in our experiments and numerical simulations. 
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I. INTRODUCTION 



^ I Information transfer through nonhnear systems in the presence of fluctuations has been 

extensively studied in the context of stochastic resonance (SR) The frontier of this 
research has shifted towards systems with spatial degrees of freedom over the past few years. 
While the efforts initially were directed towards enhancing the basic SR effect 0, recent 
work has demonstrated that noise can also sustain wave propagation . In early studies, 
Jung et al [R| showed that noise can sustain spiral waves in a caricature model of excitable 
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media. The assisting role of noise for one- and two-dimensional, nonlinear wave propagation 
was experimentally confirmed by Kad'ar et al. in a chemical medium and by Locher et 
al. m in an array of coupled electronic resonators. Noise enhanced propagation (NEP) for 
periodic signals was explored by Lindner et al. in a chain of coupled, overdamped bistable 
oscillators. The authors of Ref. furnish evidence for self-organized criticality underlying 
the creation and propagation of waves by noise in a chemical subexcitable medium. So far, 
all experiments and simulations on NEP utilized local and additive noise. To our knowledge, 
no comparative studies on the effectiveness of local vs. global, additive vs. multiplicative 
noise have been attempted. 

In this paper, we investigate experimentally the effectiveness of global noise as compared 
to local noise for the propagation of a signal in a chain of coupled bistable elements. Ex- 
panding on earlier reports the experiments are performed using a 16 x 16 array of diode 
resonators driven in the stable period-2 regime. A bias consisting of a second drive at half 
the main frequency renders one phase more stable and a phase kink can be made to propa- 
gate across the array. For an intermediate value of the coupling resistors and small bias, we 
observe propagation failure. Adding noise to the drive of each element, either from a single 
source or from individual sources - corresponding to global and local noise, respectively - 
greatly increases the chances of successful signal transmission. Digital simulations of NEP 
in a simple model of cardiac tissue compare favorably with these observations and allow us 
to investigate the effects of parameter mismatch between the elements. At last, we present 
theoretical insights borrowed from discrete kink statistics which describe the phenomena in 
reasonable agreement. 

II. EXPERIMENTAL RESULTS 

The experimental setup consists of 256 coupled diode resonators, of which one element 
is shown in Fig.^. The elements are arranged in a 16 x 16 array with periodic boundary 
conditions connected along two opposite edges. Each diode resonator works as a bistable 
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element when driven in its period-2 state. We break the phase symmetry by adding to 
the drive a second sinusoidal signal at half its frequency. This bias renders one phase more 
stable. We refer to the less stable phase as the metastable state. By inducing a phase change 
at one edge of the array a one-dimensional wavefront comprised of phase kinks will travel 
towards the detector at the other edge. The noise generators were constructed using the 
shot noise generated by a current through a pn junction diode as a source. 

It is clear that in the absence of noise, a local phase jump will lead to a "domino 
effect" only if the bias and coupling are strong enough. The energetically lower phase then 
propagates into the metastable phase in the form of a moving kink. For identical elements, 
the speed of this moving interface depends on both the coupling strength and the amplitude 
of the applied bias. If the latter two parameters are chosen low enough, kinks in discrete 
systems will fail to propagate. In the experiment, there is a third factor contributing to kink 
trapping, namely heterogeneity of the chain. In our system, the variation of key parameters 
of the diode resonators and the difference in local noise power is as high as 10 %. Arranging 
the array with periodic boundary conditions in one dimension and preserving the motion 
of wavefronts in the other dimension, we considerably reduce inhomogenities by effectively 
averaging over 16 elements. 

Figure. ^ displays the measured kink velocities in the absence of noise as a function of 
bias for an intermediate value of the coupling resistors. The velocity decreases approximately 
linearly with the bias down to a cutoff value of 0.9 units where it rapidly falls off to zero. We 
operate the array at a bias of 0.6 units for which the system is not capable of deterministic 
kink motion. 

Our experimental results are given in terms of the arrival-probability curves shown in 
Fig. 1^. The procedure for obtaining these are analogous to those described in |^ : We reset 
all resonators to be in the metastable state and then induce a phase flip at one edge of the 
array. Any successful kink propagation generates a step function at the detector, which we 
average over approximately 100 events, resulting in a smooth rise curve (solid lines). In 
order to quantify transmission degeneration by noise nucleated spurious signals we repeat 
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the same experiment without inducing a phase flip initially (dashed lines) . 

The arrival-probability curves for both local and global noise at flve different noise 
strengths are shown in Fig. ^. For local, i.e. independent from site to site, noise we 
obtain results analogous to those previously reported ||^. For a local noise background of 
less than 0.0025 mV'^/Hz (Fig. |^(a)) kinks remained trapped longer than the measurement 
time of 5000 drive cycles. At 0.0121 mV'^/Hz the signals arrive with a wide distribution of 
travel times and a slow mean (Fig. 0(b)). The arrival times become shorter with increasing 
noise strengths (Figs. |^(c) and §(d)). Finally, at 0.0625 mV'^/Hz a substantial number of 
false starts (dashed lines) corrupt the detection of the original input signal. 

We found that the qualitative behavior of the chain under the influence of global noise is 
similar, but the onset of detectable kink propagation is found at much lower noise strengths. 
We observe no kink propagation below noise levels of 0.001 mV'^/Hz. Slow and disperse 
kink motion occurs for noise levels of 0.0025 mV'^/Hz (Fig. |^(a)). Higher average kink 
speeds and less fluctuations are recorded for levels of 0.0121 and 0.0256 mV'^/Hz (Fig. ^(b) 
and 0(c)), while the signal is severely corrupted for values > 0.04 mV'^/Hz (Figs. |^(d) and 
1(e)). 

Fig. ^ compares the velocities of the propagating wavefront as a function of the noise 
strength for both global and local noise. The velocity of the propagating wavefront shows 
an approximate linear increase with increasing noise strength in both cases. Note that 
the velocity is simply the inverse of the measured arrival times -multiplied by number of 
sites. It is hence well-deflned only for low noise levels; in the case of substantial nucleation of 
additional thermal kinks-antikink pairs this calculated "velocity" should be interpreted with 
caution. Therefore, the data points in Figs. | and |^ which correspond to signiflcant noise 
corruption should be considered as outliers. For the coupling resistors used the velocity for 
the global noise case is about 15 % greater than that seen for local noise. 

Besides the earlier onset in the case of global noise, there are also notable differences in 
the mechanism that leads to spurious signals. Clearly, for identical elements - unlike spatially 
uncorrelated noise - global noise cannot induce spurious kinks, which then compete with the 
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deterministic kink due to the signal. The only way of generating a "false alarm" at the 
detector would be to phase flip the entire chain, thus requiring a single large fluctuation. 
This picture is not fully correct if there are mismatches between the elements but illustrates 
the dominant mechanism of creating the signal-masking noise background. Independent 
noise sources, however, easily spawn spurious kinks but rarely ever cause a global phase flip 
of the array. We postulate that the mechanism of noise sustained propagation per se is not 
very different for local and global noise, presumed the kink width is small (i.e. involving 
a few elements only). The reason for this conjecture is the local nature of the stochastic 
escape processes that provide for the average effective kink displacement. As long as the 
spatial correlation length of the noise is not substantially smaller than the kink width, the 
noise induced "propulsion " of the kink will be similar for global and local noise. Note that 
the findings of Kadar et al. 0] who briefly addressed the issue of noise correlation lengths, 
are consistent with this conjecture. 



III. SIMULATIONS OF A MODEL OF CARDIAC TISSUE 

Propagation failure of signals due to discreteness of the supporting medium has pre- 
viously been observed in theoretical and experimental studies of cardiac tissue. In 
particular. Keener introduced a modified cable theory, which incorporates the discretizing 
effects of the so-called gap junctions 0. Gap junctions, characterized by the (relatively 
high) intercellular resistance r^, provide the electrical coupling between cardiac cells. Math- 
ematically, the propagation of action potential along cardiac cells is described by various 
cable theories, which are analogous to wave propagation in one-dimensional conductors (ca- 
bles). Formally, continuous and discrete models describe the wave propagation by either a 
partial differential equation or - in the latter case - via coupled ordinary differential equa- 
tions. Continuous cable theory either ignores the effects of the gap junctions or replaces the 
cytoplasmic resistance with an effective resistance; in either case the electrical resistance is 
assumed to be spatially homogeneous. Here, we focus on the opposite assumption that gap 
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junctional resistance is much more important than cytoplasmic resistance. We thus neglect 
the dynamics within a cell and assume that the propagation of the action potential is domi- 
nated by the delay caused by the gap junctions 0. Within the context of this discrete cable 
theory, we can write the current balance as 

CmS'^ = — {(pn+l - 20„ + + SIm{(l)n) (l) 

where 0„ is the transmembrane potential for the nth cell, S the surface area of cell membrane, 
and Cm is the membrane capacitance per unit area of membrane. specifies the inward 
ionic currents per unit area of membrane and is generally postulated to be a function of 
(pn having three zeros. For simplicity, we choose a simple cubic polynomial SI„i{(t>) = 
12v^0(l -<p){<p- 0.5) + 0.5 [|. Note that though similar in appearance, Eqn. (|l]) is not 
simply a discretization of its continuous analog 

(where L is the size of the cardiac cell), but stands in its own right as a spatially discrete 
nonlinear wave equation. The most important observation is that propagation can fail in 
model (|I]) if is sufficiently large, but increasing resistances in Eq. (Q) can never lead to 
propagation failure. Note that depends on the excitability of the tissue. 

Fig. ^ shows a plot of the numerically determined speed (solid line) of propagation for 
model (P as a function of the coupling strength d = l/r^ as well as the analytically obtained 
kink speed c = ^'^"^ \fd (dashed line) for Eq. (Q). The reader is referred to Ref. for an 
explanation of cq and Rm- It is evident that propagation is impossible for Vg larger than a 
certain critical value r*, which turns out to be a monotonic increasing function of excitability 

i- 

We have performed digital simulations of the stochastic modification of ([^) 

^ = e (0„+i - 20„ + </)„„i) + 5/™(0„)(l + 6/(i)) + U{t) (3) 

using the Euler-Maruyama algorithm with a time-step of dt = 0.05 and coupling strength 
e = 0.07 (40 elements), ^a(^) and ^Af (^) are additive and multiphcative Gaussian white noise, 
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bandlimited in practice by the Nyquist frequency In = We quantify the noise by its 
dimensionless variance = 2DfN, where 2D is the height of the one-sided noise spectrum. 

Here, we only consider the case of purely additive noise, ^Af(^) = 0. By following an 
analogous procedure as in the experiment we obtain the probabilities for successful signal 
transmission as illustrated in Fig. ^. As in the experiment, global noise provides for kink 
propagation at much lower values of o"^ than local noise does. Fig. ^ compares the velocities 
of the noise propulsed kink as a function of for both global and local noise. The velocity 
of the propagating wavefront shows an approximately parabolic dependence on the noise 
power in both cases. For the coupling strength employed global noise leads to speeds about 
15 % greater than that observed for local noise. 



IV. DETECTION CRITERIA 

The decision whether a kink arriving at the last element corresponds to a signal injected 
at the first site constitutes a simple binary hypothesis-testing problem |T0|. We assign the 



null hypothesis Hq to "no signal injected" and the opposite for the alternative hypothesis 
Hi. Denote the according decisions Di as the judgment that hypothesis Hi was in effect. 
Clearly, there are two possible errors: A so-called Type I error occurs when making the 
decision that Hi was in effect while the contrary is true. Borrowing notation from radar 
detection, we refer to this as the probability of false alarm Pf = P{Di\Hq). On the other 
hand, if Hi was in effect to generate the data, and we decide Dq, then we have committed 
a Type H error, which in radar is referred to as a missed detection We do that with 
some probability P{Dq\Hi) which is related to the probability of detection P^ = P{Di\Hi) 
in an inverse fashion: Pa = 1 — P{Dq\Hi). 

In our experiment, there is no uniquely defined, objectively "best" noise level without first 
defining a decision strategy. Two suitable approaches are (i) the N eyman-Pearson strategy, 
in which the probability of detection Pd is maximized while specifying an upper bound 
for the false alarm probability Pf, and (ii) Bayes' rule which assigns costs to the various 
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outcomes of the decision process, such as correctly detecting a signal or being "deceived" 
by a spurious kink. The optimum noise level in the latter case would be the one, which 
minimizes the total average cost. 

In communication systems one is usually interested in the total probability of error Pg = 
P{Di\Hq) + P{Dq\Hi) = Pf + 1 — Pd- Though in the experimental setup there is no 
inherent time-scale, i.e. the decision when to reset the chain is rather arbitrary, in digital 
communication applications we would expect information bits to be sent at a constant rate. 
Hence, we choose a reasonable time interval at the end of which we measure the probabilities 
of false alarm Pf and missed detection 1 — as a function of noise power. Then, Pf is 
simply the value of the dashed line in Fig. ^ at time = 300, and P^ is the corresponding value 
for the solid line. For local and global noise, the total probability of error Pg is displayed in 
Fig. for very low and rather high values of the noise power, Pg is almost one. However, 
there exists an optimal noise strength for which both Pf and 1 — Pa nearly vanish, resulting 
in a sharp minimum of Pg. 

V. THEORY 

The results of Sects. || - |V| touch upon two central properties of kink statistics in a 
discrete bistable chain, namely the Brownian motion of an individual kink and the nucleation 
of kink-antikink pairs in the presence of an external static bias (or d.c. forcing term). In 
the absence of a precise model for the 2-state potential that describes the phase shifts 
in the diode resonators of Sec. H or drives the transmembrane potential of Sec. |T|, we 
speculate that the in situ bistability of our arrays can be rendered satisfactorily in terms 
of a Double Quadratic (DQ) potential, that is by two parabolas displaced by a distance 2a 
(see Fig. ||a). A discrete DQ model is likely to capture, at least qualitatively, the essential 
features of the array dynamics investigated above, while affording substantial simplifications 
in its analytical treatment. The analysis of this Section can be carried over, with more 
mathematical effort, to the 0^ model of Sec. as well. 
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The DQ model has been studied both in the continuum [|TT| and in the discrete case [|T3 . 
In dimensionless units the DQ Hamiltonian reads 

^ = { Y + - + (0. - + ^(101 - 1)^} , (4) 

with Hq = ma'^/l. Each 0„ can be regarded as the displacement (in units of 2a) of the n-th 
chain site with mass m, Cq and ujq represent respectively the limiting speed and frequency 
of the phonon modes propagating along the chain and / denotes the chain lattice constant ( 
for instance, in Eqn. (|l|) I was set to one). The ratio Cq/P, which quantifies the effectiveness 
of the coupling between two adjacent bistable units, is the couphng constant of our model. 
The importance of the discreteness effects is measured by the discreteness parameter 

namely the ratio of the kink length d to the chain spacing I. 

A. The continuum limit 

In the continuum (or displacive) limit 7 oo the Hamiltonian (^ can be expressed as 
the line integral of the Hamiltonian density 

^M = y + c^f + (6) 

2 

where the string field (p{x,t) is defined as lim/^o 0a:/«(^) ^[4>] = "^(101 "~ 1)^- The 
statistical mechanics of the continuum DQ model can be worked out analytically in great 



detail [|rT|. In particular, we know that the kink (0+) and the antikink solutions ((/)_) 

0±(x,t) = ±sgn[x - X{t)][l - expV^-"'/=o] (7) 

can be regarded as relativistic quasi-particles with size d = Cq/uq, mass Mq = Eq/cq = 1/d 
(or rest energy Eq = uqCo) and center of mass X(t) = Xq + ut. 

At low temperatures, kT <^ Eq, any string configuration can be represented as a linear 
superposition of randomly distributed kinks and antikinks floating on a phonon bath. A DQ 



string in equilibrium at temperature T and with boundary conditions 0(— oo, t) = 0(+oo, t) 
bears naturally a dilute gas of thermal kink-antikink pairs with density 

The qualification thermal underscores the fact that no pairs per unit of length (with nod <^ 1) 
are being generated by thermal fluctuations alone, irrespective of any geometric constraint 
at the boundaries (see discussion of Fig. ^ in Sec. ||). 



B. Kink Brownian motion 

The 0±(x,t) solutions (0) tend to travel with arbitrary constant speed u < Cq, unless 
perturbed by the coupling to a heat bath or by an external field of force. The simplest heat- 
bath model was obtained 0] by adding a viscous term —a(j)t and a zero-mean, Gaussian 
local noise source ({x, t) to the string equation of motion corresponding to the Hamiltonian 
density (H), that is 

(i>tt - cl(t),x - u;osgn[0](|0| - 1) = + F + C{x,t). (9) 

Note that all our experiments and simulations have been carried out in the overdamped limit, 
a ^ tuo; and in the presence of an additional sub-threshold force F with F < uJq, also incor- 
porated in Eq. (|^). Thermalization is imposed here by choosing the noise autocorrelation 
function (C(x, t)C{x\ t')) = 2akT5{t - t')6{x - x'). 

A single kink (antikink) subjected to thermal fluctuations undergoes driven Brownian 
motion with Langevin equation 

2F 

^ = T-rr+^W, (10) 

where rj{t) is a zero-mean- valued Gaussian noise with strength D = kT/aM^ and autocor- 
relation function {7]{t)r){0)) = 2D6{t). As apparent from Eq. (|10[), the external bias pulls 
(j)± in opposite directions with average speed ±uf and up = 2F/aMo. 
If the local fluctuations are spatially correlated, say 
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{C{x,t)C{x',t')) = 2afcT5(t-t')[e"'"""''/V2A], (11) 



the noise strength D changes into [14 



... Dd f \ \ 

= 1 1 + ^ 1 ■ (12 

As speculated in Sec. 0, for noise correlation length A smaller than the kink size d possible 
spatial dishomogeneities become negligible, i.e. -D(A) ~ D for A ^ d [1^. The global noise 
regime simulated both experimentally in Sec. ^ and numerically in Sec. |ITT| corresponds to 
the limit A — oo of the source ({x, t) rescaled by the normalization factor v^2A; the Langevin 
equation ( |I0| ) still applies, but the relevant noise strength is now limA-+oo 2AD(A) = AD. This 
accounts for the observation that global noise sustains kink propagation more effectively than 
local noise. Note that the enhancement factor of 4, more exactly 4a, is nothing but twice 
the distance of the DQ potential minima (in dimensionless units). 

Another important property of global noise is that it cannot trigger the nucleation of a 
kink-antikink pair and, therefore, minimizes the chances of a "false alarm" (see Fig. |^). For 
this to occur it would be necessary that a spatial deformation of a stable string configuration 
(vacuum state) be generated large enough for the external bias to succeed in making it grow 
indefinitely. Such a 2-body nucleation mechanism would require a /oca/ breach of the (f) — </> 
symmetry of the DQ equation (|^), which can be best afforded in the presence of uncorrelated 
in situ fluctuations ||T6| . 

The nucleation rate, namely the number of kink-antikink pairs generated per unit of time 
and unit of length, can be easily computed by combining the nucleation theory of Ref. 



with the analytical results of Ref. |TI|] for the DQ theory. For values of the string parameters 
relevant to Sees. || - that is for kT and Fd -C Eq the stationary DQ nucleation rate can 
be approximated to [0 

ri(T) = = 2u,nl{n (13) 

if Fd < kT, or T2{T) = iy^ri(T), if kT <^ Fd <^ Eq 0. For an overdamped string, 
a ^ ti^o the time constant t(T) amounts to the kink (antikink) lifetime prior to a destructive 
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collision with an antikink (kink). Both estimates for the DQ nucleation rate clearly show 
that spontaneous nucleation of thermal pairs may degrade appreciably local-noise sustained 
propagation of injected (or geometric) kinks only for thermal energy fluctuations of the order 
of the kink rest energy. 



C. The Peierls-Nabarro potential 

Let us go back now to the case of a discrete DQ chain. Discreteness (with parameter 7) 
affects the kink dynamics on two accounts: 

(i) The proflle of a static kink (antikink) 0±(a;, 0) is deformed into 



0g„ = ±sgn[n - iV] [1 - Z.z/l-^l], (14) 

with = 2y^/l + u, N = m + 1/2, m = 0, ±1, ±2, . . ., and z/ = [^/TT4^ -l]/[VlT4^ + 
1]. To make contact with the displacive solution 0±(x, 0) one must replace nl with x, Nl 
with Xq and take the continuum limit 7 ^ 00 (so that z/ ~ 1 — I/7). Note that the spatial 
extension of the discrete kink solutions (j)±]n increases monotonically with 7. As 7 decreases 
below unity, approaches a step function (order-disorder limit); 

(ii) (j)±]n is centered midway between two chain sites due to the conflning action of an 
effective [or Peierls-Nabarro (PN)] potential [|13|. The PN potential describes the spatial 

(s) 

modulation of the (f)±'n rest energy as its center of mass is moved across one chain unit cell, 
say from ml up to (m + 1)/. 

As a result, according to the Langevin equation approach of Sec. |VB| , the (j)±]n center of 
mass X{t) diffuses on a periodic, piece-wise harmonic potential with constant I and angular 
frequency upn, that is [1^ 



aX = -uj'pj^lX - /(Int[X//] - 1/2)] t 2F/Mq + ar/(t), (15) 

where cup^y ~ (1 + v)uJq and Int[X//] denotes the integer part of X in units of /. Note that 
^PN ^0 and cupTv — > \/2ujq in the highly discrete and continuum limit, respectively. The 
energy barriers of the PN potential are thus (almost) quadratic in /. 
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The one-dimensional Langevin equation ([Tq ) has been studied in great detail by Risken 



and coworkers ||T^. In the noiseless limit, ri{t) = 0, the process X{t) is to be found either 
in a locked state with (X) = 0, for AF/Mq < ujpjy , or in a running state with (X) ~ up, 
for AF/Mq > ujpj\f. This is indeed the depinning (or locked-to-running) transition described 
in Fig. ^ At finite temperature the stationary velocity (X) = u{T) can be cast in the form 
following, 

^ = 1 (16) 

where 6 = 2FI / kT and the quantities A and B can be computed numerically with minimum 
effort |jl8|. The ratio u{T)/up is the rescaled mobility; it crosses from (locked state) 
over to 1 (running state) continuously in a relatively narrow neighborhood of the threshold 
value Fth = moreover, u{T)/up increases monotonically with T at fixed bias. 

Such a temperature dependence of the kink mobility explains the sequences of rise curves in 
Figs. ^ and |^, where kink propagation seems to speed up on raising the noise level. 

VI. SUMMARY 

In conclusion, the present analysis confirms our speculation that the apparent SR behav- 
ior of the efficiency of noise-sustained transmission of kink-like signals along a bistable chain 
results from two competing mechanisms, both controlled by noise: The driven diffusion dy- 
namics of stable noninteracting kinks, which increases exponentially with the temperature in 
the vicinity of the depinning transition (propulsion mechanism); The detection of spurious 
signals, as thermal kink-antikink pairs nucleate with exponentially increasing rates, thus 
corrupting the propagated signal (garbling mechanism). 

If the spatial distribution of the noise was constrained to a small neighborhood around the 
kink and zero along the rest of the chain, fast and efficient noise supported signal transmission 
without false alarms would be realizable. This seemingly artificially constructed scenario can 
be achieved naturally by considering the case of purely multiplicative noise [IT3. A detailed 
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study of noise sustained propagation in the presence of multiplicative fluctuations is beyond 
the scope of this work. 

During the preparation of this manuscript the authors learned about recent results on 



propagation failure in the context of cell differentiation |^. Utilizing a highly simplified 
model composed of coupled bistable elements, the authors furnish evidence for the discrete 
nature of chemical signaling waves propagating through a chain of cells. We speculate that 
fluctuations, inherent in biological systems, might play a significant role in the details of cell 
differentiation processes. 

We acknowledge the Office of Naval Research for financial support. ML, NC and EH 
warmly thank D. Cigna for very significant contributions to the experimental setup. 
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FIG. 1. Experimental circuit arrangement for any node (i,j). We either add noise locaUy, as 
shown, or globally by adding one noise source to the main drive. Coupling is provided by the 
resistors Rc- 

FIG. 2. Kink velocity as a function of the bias for an intermediate value of the coupling resistor 
Rc- The kink speed shows an an approximate linear decrease with bias down to a cut off value of 
approximately 0.9 units. 

FIG. 3. Arrival-probablity of a signal at the end of the chain with (solid line) and without 
(dashed curves) induced kink at the beginning, for five local and global noise levels increasing from 
top to bottom. 

FIG. 4. Experimentally observed average kink velocities as a function of noise strength for both 
global and local noise. Note that the last (two) data point (s) for local (global) noise are corrupted 
by noise and should be interpreted as velocity only with caution. 
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FIG. 5. Kink speed as a function of 
coupling strength d = — for the con- 
tinuous model Eqn. (P) (dashed line) 
and the discrete model Eqn. (Q) (solid 
line). Note that propagation in the dis- 
crete model is impossible for rg > r* 
{d < d* = p^). For this simulation, CmS = 1 
and SImi(p) = 12\/3(/)(l - </))(</) - 0.5) + 0.5 

FIG. 6. Arrival probabilities for the discrete model Eqn. (|^) with (solid line) and without 
(dashed curves) induced kink at site 1, for five local and global noise levels increasing from top to 
bottom. The propagation distance spans 40 elements. Note that we employ dimensionless units 
for the noise power a^. 

FIG. 7. Numerically measured average kink velocities as a function of noise strength (dimen- 
sionless)for both global and local noise. It is insightful to compare these values with Fig. ^. As in 
Fig ^, the last data point includes a large fraction of spurious kinks which results in an artificially 
high value for the velocity. 

FIG. 8. Total probability of error Pg as a function of noise variance o"^ (dimensionless) for global 
(solid line) and local noise (dashed line). For a range of optimal noise strengths Pg virtually 
vanishes. The qualitative behavior is robust to variations in the measurement time, which here is 
taken to be 300 time units. 
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FIG. 9. The Double Quartic (DQ) model 
for LjQ = cq = a = 1 and discreteness parame- 
ter 7 = 1. (a) The DQ potential V[(j)] of Eqs. 
(D and (|); (b) The PN potential V{x,F) 
of Eq. (ll|) and ^ for F/Fth = 0.08; (c) 
Kink stationary velocity u{T) versus bias in- 
tensity F for Wq/ZcT = 20 (curve 1) and 50 
(curve 2). At T = 0+ the limiting curve 
(dashed) is given hy u/up = for F < Fth 
and u/uF = {(F/Fth) ln[F/(F - Ft^)]}-^ for 
F > Fth- All quantities plotted are dimen- 
sionless. 
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